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Influence of shower fluctuations on the shape of lateral distribution of electrons in EAS of fixed size measured 
by scintillation counters is analyzed in framework of scaling formalism. Correction factors for the mean square 
radius of electrons are calculated for the experimental conditions of KASCADE array. Possible improvement of 
the primary mass discrimination by analysis of lateral distribution of EAS electrons is discussed in detail. 



1. INTRODUCTION 

The determination of the primary cosmic ray 
chemical composition from extensive air showers 
(EAS) observations is an open problem. Vari- 
ous techniques based on different EAS character- 
istics measured by different experiments including 
multi-component methods do not exhibit consis- 
tent results in estimation neither primary mass in 
the case of individual showers nor the mean mass 
composition at a certain energy. 

One of the key EAS quantities necessary for ba- 
sic shower parameters reconstruction is the lateral 
distribution of charged particles at fixed observa- 
tion depth. The exact form of lateral distribution 
function (LDF) is still uncertain. The majority 
of analytical parameterizations of LDF of differ- 
ent EAS components is traditionally based on the 
well known Nishimura-Kamata-Greizen (NKG) 
function pQ: 
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Here p(r;i?o,s) is the particle density at radial 
distance r from the core position in shower with 
primary energy Eq and the age parameter s, 
N(Eq,s) - total number of particles at the ob- 
servation depth, Rq - shower scale radius, which 
does not depend on primary particle type and 
energy (originally - the Moliere unit). Various 
modifications of NKG form, such as introducing 



an additional fixed or age-dependent scale coeffi- 
cient or a local age parameter s(r) and also gen- 
eralizations of the function by using third power- 
law term were suggested. 

A different theoretically motivated approach 
(scaling formalism) was proposed in our pa- 
pers 03]: 
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where t is the observation depth. The scaling 
function F(X) can be described as follows: 



F(X) = CX~ a (1 + X) 
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For electron densities we find C — 0.28, a — 1.2, 
j3 — 4.53, 8 — 0.6, Rq = R m . s . - root mean square 
radius of electrons: 
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According to our calculations, the scaling for- 
malism allows to reproduce electron LDF with 
10% uncertainty for E„ = (10 14 - 10 20 ) eV, 
t = (600 - 1030) g/cm 2 , X = (0.05 - 25). The 
last condition corresponds to the radial distance 
range from r ~ (5 — 10) m tor ~ (2.5 — 4) km de- 
pending on the shower age. This limitation makes 
scaling approach inadequate for shower size and 
core position estimation, but accurate enough for 
description of the shape of LDF measured by 
ground-based shower arrays far from the core. 
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Figure 1. Primary energy distributions in four selected bins N e , simulated according to KASCADE shower 
classification procedure [6J (distribution are renormalized to equal height of maximum for convenience of 
comparison). 



The method for mean primary mass deduction 
based on scaling formalism was developed in |4l5j . 
The advantage of this method is its applicability 
to the experimental data of both compact and 
giant air shower arrays in wide primary energy 
range and also relatively weak sensitivity of the 
conclusions to variations of basic parameters of 
hadronic interaction model implemented in calcu- 
lations. Unfortunately, the shape of charged par- 
ticle LDF measured experimentally is affected by 
the experimental method of shower classification. 
This effect can also be described in the framework 
of scaling formalism, but the variation of the root 
mean square radius compared with the data ob- 
tained theoretically for the fixed primary energy 
should be taken into account in case of relatively 
low primary energies, when shower selection is 
made by the total number of electrons, e.g. for 
example for KASCADE and Moscow State Uni- 
versity air shower arrays. 

In this paper we examine thoroughly the influ- 
ence of shower fluctuations on the shape of lat- 
eral distribution of electrons in EAS of fixed size 
under the experimental conditions of KASCADE 
array in order to improve the reliability of mean 
primary mass deduction from the experimentally 
measured LDFs. 



2. CALCULATION METHODS 

We made simulations of extensive air showers 
initiated by protons and iron nuclei of vertical 
incidence assuming power-law differential energy 
spectrum of primaries with exponent ol\ = 2.62, 
a2 = 3.02 and also with sharp knee from ai 
to a 2 at Eq = 10 6 - 5 GeV. We used the semi- 
analytical code [5] with full Monte-Carlo treat- 
ment of hadronic part of cascade based on quark- 
gluon string model and analytical expressions of 
pure electromagnetic subshowers keeping all the 
basic sources of fluctuations. The fluctuations of 
different EAS components calculated by our code 
are in good agreement with CORSIKA/QGSjet 
results. 

According to |6J we simulated shower classifi- 
cation procedure used at KASCADE array and 
evaluated lateral distributions and root mean 
square radiuses of electrons in eight bins of shower 
size. The number of showers in each bin amounts 
from ~ 5000 for lower energies (lg N e = 3.9 — 4.3) 
to - 1000 for higher energies (lg N e = 6.7 - 7.1). 

3. RESULTS 

The primary energy distributions in four from 
eight bins of shower size is shown in fig. [1] One 
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Figure 2. Correction factors K — (i^s/^ms) f° r 
proton initiated vertical EAS at sea level assum- 
ing different primary energy spectrum exponents: 
ct\ = 2.62 (triangles with dotted approximation 
curve), o;2 = 3.02 (squares with dashed curve), 
spectrum with the knee at Eq = 1 6 5 GeV (solid 
circles with solid curve). See text for details. 



can see that the energies largely overlap in differ- 
ent bins though the selected bins are not neigh- 
boring. 

We evaluated the correction factors defined as 
the ratio of root mean square radius calculated for 
certain shower size bin to that for corresponding 
average primary energy: K = (R^/R^), These 
correction factors for vertical proton initiated 
showers at sea level calculated with above men- 
tioned assumptions about primary energy spec- 
trum are shown in fig. [51 It is clear, that values 
of K approach to 1 with energy as shower fluctu- 
ations decrease. The same effect takes place for 
a heavier primary nuclear or for smaller shower 
size bins. At the same time the correction, 
which should be made for an adequate compar- 
ison of theoretical and experimentally estimated 
root mean square radiuses is essential for all con- 
sidered shower size bins. 

Assuming the validity of scaling approach we 
made one-parametric fitting of experimental lat- 
eral distribution of electrons obtained by KAS- 



Figure 3. The superposed radial scale factors Rq 
obtained from experimental data of KASCADE 
under different assumptions about scaling func- 
tion in comparison with root mean square radius 
calculated for vertical proton and iron initiated 
showers at sea level. 



CADE array [6] using formula ([2]) with three dif- 
ferent scaling functions: 

1) theoretically proved function ([3]) with mean 
square radius of electrons as radial scale parame- 
ter; 

2) modified NKG-function [5] with fixed shower 
age parameter (s=1.65) and variable Rq; 

3) polinomial function: 

F(X) = Cexp j^^mX)*!, (5) 

with n = 4 and parameters at being fixed for 
all bins independently from Rq. Some discrimi- 
nation of experimental data at small radial dis- 
tances was done in order to eliminate points with 
r < 0.05i? ms . 

All the fitting functions give satisfactory overall 
fit of experimental data of KASCADE with resid- 
uals not exceeding 10% for polynomial and mod- 
ified NKG functions and 15% for function ([3|). 
Though polynomial and NKG functions give con- 
siderably better accuracy in considered radial dis- 
tance range, they both lead to incorrect predic- 
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Table 1 

Radial scale factors (Rq ±SRq) obtained by fitting of experimental LDF [5] by different scaling functions. 





Function © 


Modified NKG (s=1.65) 


Polinomial © 


3.9-4.3 


146.8 ±2.5 


29.79 ± 3.0 • 10" 1 


174.4 ± 1.6 


4.3-4.7 


134.3 ±2.2 


26.93 ±2.1 ■ 10" 1 


156.6 ± 1.2 


4.7-5.1 


125.0 ± 1.5 


25.17 ± 1.9 ■ 10" 1 


146.0 ± 1.2 


5.1 - 5.5 


122.6 ± 1.3 


23.92 ± 1.3 • 10- 1 


138.6 ±0.9 


5.5 - 5.9 


122.8 ± 1.6 


23.65 ± 1.1 ■ 10" 1 


137.8 ±0.5 


5.9 - 6.3 


122.6 ±2.1 


24.03 ± 1.6 ■ 10" 1 


139.8 ±0.9 


6.3 - 6.7 


125.4 ±2.3 


24.70 ±2.1 ■ 10" 1 


141.0 ± 1.2 


6.7-7.1 


130.6 ± 1.8 


25.55 ± 3.5 • 10" 1 


145.1 ±2.1 



tions for extremely large core distances, while the- 
oretically motivated function ([3]) remains realistic 
up to r ~ 25i? ms . 

The values of radial scale parameters i?o ob- 
tained by fitting the experimental data are sum- 
marized in Table [T] It is not surprising, that dif- 
ferent fitting functions correspond to significantly 
different values of Rq. An additional bias in Rq 
can be related with the insufficiency of radial dis- 
tance range well covered by the array or some 
other systematic errors in data processing. So 
it is worth to compare the rate of change of ra- 
dial scale factors with energy <9i? / d log Eq which 
obviously reflects the rate of change of mean pri- 
mary mass. In fig. Rvalues of i?o superposed with 
each other by the appropriate factors are shown in 
comparison with root mean square radius calcu- 
lated for vertical proton and iron initiated show- 
ers at sea level taking into account the correction 
factor K for primary spectrum with the knee. As 
it is seen from the figure the rate of change of Rq 
obtained using different functions is consistent. 

4. CONCLUSIONS 

1. Correction by a factor K = (R^' s /Rf ns ) 
should be made when comparison of lateral dis- 
tributions of electrons measured by ground-based 
experimental arrays to LDFs calculated theoreti- 
cally for fixed primary energy is carried out. 

2. Absolute values of radial scale factor R con- 
tain systematic errors and could be biased de- 
pending on the form of lateral distribution func- 



tion chosen for experimental data processing and 
final fitting. However, if one concerns the rate 
of change of R Q with primary energy, which is a 
good measure for primary mass composition vari- 
ation, then different scaling functions used in the 
framework of scaling formalism do not contradict 
each other. 

3. Basic analysis using different assumptions 
about the form of scaling LDF leads to consistent 
model insensitive conclusion that average primary 
particle mass above the knee increases with en- 
ergy, that is in reasonable agreement with the 
large number of experimental observations (see 
e.g. [7]) and also with recent results of the anoma- 
lous diffusion model [5]. 
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